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ABSTRACT 

We present continuous, monochromatic star formation rate (SFR) indicators over the mid-infrared 
wavelength range of 6 — 70 ^m. We use a sample of 58 star forming galaxies (SFGs) in the Spitzer- 
SDSS-GALEX Spectroscopic Survey (SSGSS) at z < 0.2, for which there is a rich suite of multi¬ 
wavelength photometry and spectroscopy from the ultraviolet through to the infrared. The data from 
the Spitzer infrared spectrograph (IRS) of these galaxies, which spans 5 — 40 pm, is anchored to their 
photometric counterparts. The spectral region between 40 — 70 pm is interpolated using dust model 
fits to the IRS spectrum and Spitzer 70 and 160 pm photometry. Since there are no sharp spectral 
features in this region, we expect these interpolations to be robust. This spectral range is calibrated as 
a SFR diagnostic using several reference SFR indicators to mitigate potential bias. Our band-specific 
continuous SFR indicators are found to be consistent with monochromatic calibrations in the local 
universe, as derived from Spitzer, WISE, and Herschel photometry. Our local composite template 
and continuous SFR diagnostics are made available for public use through the NASA/IPAC Infrared 
Science Archive (IRSA) and have typical dispersions of 30% or less. We discuss the validity and 
range of applicability for our SFR indicators in the context of unveiling the formation and evolution 
of galaxies. Additionally, in the era of the James Webb Space Telescope this will become a flexible 
tool, applicable to any SFG up to z ~ 3. 

Subject headings: galaxies: star formation — infrared: galaxies — stars: formation 


1. INTRODUCTION 

Star formation is a fundamental parameter of galaxies 
that describes how galaxies evolve, when used in conjunc¬ 
tion with mass. As the process of star formation depletes 
a galaxy of its gas, it must be continuously replenished 
by infall from the intergalactic medium to be supported 
for an extended time. When massive stars die, they 
enrich the surrounding interstellar medium with heavy 
metals, thus altering a galaxy’s chemical composition. 
Therefore, accurately tracing star formation though cos¬ 
mic time gives key constraints on how galaxies are able 
to form and evolve (e.g., Tinsley 1968; Somerville et al. 
2012; Madau & Dickinson 2014, and references therein). 

For these reasons, great efforts have been made to cali¬ 
brate a wide range of the electromagnetic spectrum that 
can be linked to processes involved with recent star for¬ 
mation (see review by Kennicutt & Evans 2012). In par¬ 
ticular, infrared (IR) wavelength calibrations are prov¬ 
ing to be critical to understanding galaxies in the early 
Universe. Deep IR surveys with the Spitzer and Her¬ 
schel Space Telescopes have revealed that the majority 
of star formation that occurs at redshift z ^ 1 — 3 is en¬ 
shrouded by dust (e.g., Murphy et al. 2011; Elbaz et al. 
2011), making it very difficult to measure accurate SFRs 
at optical wavelengths. In addition, IR-bright galaxies 
(L > xIO^^Lq) are much more prevalent during that 
time than today (e.g.. Chary & Elbaz 2001; Le Floc’h 
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et al. 2005; Magnelli et al. 2009; Murphy et al. 2011; 
Elbaz et al. 2011; Lutz 2014, and references therein). 
Furthermore, observations suggest that ^85% of today’s 
stars were formed at redshift 0 < z < 2.5 (Marchesini et 
al. 2009; Muzzin et al. 2013; Tomczak et al. 2014). To¬ 
gether these results have renewed interest in monochro¬ 
matic (i.e., single-band) mid-IR (MIR) star formation 
rate (SFR) indicators, as distant galaxies can easily be 
observed in the MIR. 

Dust emission in the MIR is more closely related to 
star formation than longer IR wavelengths, where heat¬ 
ing by low mass (i.e., long-living) stars becomes impor¬ 
tant, which has allowed several wavelength bands in the 
MIR to be well calibrated locally as SFR diagnostics (Zhu 
et al. 2008; Rieke et al. 2009; Calzetti et al. 2010). How¬ 
ever, difficulties arise in utilizing local calibrations be¬ 
cause the regions of rest-frame wavelengths probed by a 
given band will vary with redshift. As a reference, the 
Spitzer 24 pm and the Herschel 70 pm bands target the 
rest-frame 8 pm and 23 ^m emission, respectively, for a 
galaxy at redshift z = 2. Correcting for this effect is most 
commonly achieved through A:-corrections which depend 
heavily on the assumed galaxy spectral energy distribu¬ 
tion (SED) template. Such templates (e.g.. Chary & El¬ 
baz 2001; Dale & Helou 2002; Polletta et al. 2007; Rieke 
et al. 2009; Brown et al. 2014) typically require many 
photometric bands for accurate matching. Therefore, in 
order to fully utilize current and future deep IR imag¬ 
ing surveys for a greater understanding the formation 
and evolution of galaxies without a reliance on extensive 
multi-band imaging, continuous single-band SFR indica¬ 
tors will be imperative. 

In the near future, the Mid-Infrared Instrument (MIRI; 
5-28 pm) on board the James Webb Space Telescope 
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(JWST) will expand our ability to probe galaxies in 
the MIR, detecting down to the regime of normal star¬ 
forming disk galaxies (L < 3 x out to z = 3 and 

representing an order of magnitude improvement in sen¬ 
sitivity over Spitzer bands of similar wavelength^. Thus, 
current and future cosmological surveys are highlighting 
the need for continuous monochromatic SFR indicators 
that cover, without breaks, the MIR wavelength range 
of 6 — 70 ^m. This will provide a flexible tool that can 
be applied to any galaxy up to redshift 2 ~ 3. With 
the release of the Wide-held Infrared Survey Explorer 
(WISE) All-Sky Survey (Wright et al. 2010), times are 
ripe for consolidating all these data into a coherent pic¬ 
ture. In this study, we use GALEX, Sloan Digital Sky 
Survey (SDSS), WISE, and Spitzer data of a large sam¬ 
ple of local galaxies to perform the calibration of SFR(A) 
in the 6 — 70 /rm range. 

Throughout this work we adopt the WMAP Hve-year 
cosmological parameters, Hq = 70.5 km/s/Mpc, Dm = 
0.27, Da = 0.73 (Komatsu et al. 2009). We assume a 
Kroupa (2001) initial mass function (IMF) for all SFR 
calibrations. The infrared luminosity, Air, of a galaxy 
refers to the integrated luminosity over the region from 
8 to 1000 pm (Air = /g L^dv). 

2 . DATA 

2.1. The SSGSS Sample 

The Spitzer-SDSS-GALEX Spectroscopic Survey (SS¬ 
GSS) is a sample of 101 galaxies located within the 
Spitzer Wide-Area Infrared Extragalactic (SWIRE) Sur- 
vey/Lockman Hole area at 0.03 < 2 ; < 0.22 (Treyer et 
al. 2010; O’Dowd et al. 2011). These galaxies repre¬ 
sent a subset of the 912 galaxies within the Johnson et 
al. (2006) sample, which has extensive multi-wavelength 
coverage from Spitzer^ SDSS, and GALEX, and for which 
Spitzer infrared spectrograph (IRS) measurements have 
also been obtained. The UV data is from pipeline- 
processed GALEX observations of this regions with aver¬ 
age exposures of ~1.5 ks. The optical photometry comes 
from the seventh data release of the SDSS main galaxy 
sample (DR7; Abazajian et al. 2009). The optical spec¬ 
troscopic measurements are from the Max Planck In¬ 
stitute for Astrophysics and Johns Hopkins University 
(MPA/JHU) group®, which is based on the method pre¬ 
sented in Tremonti et al. (2004). The infrared photom¬ 
etry comes from the SWIRE survey observations (Lons¬ 
dale et al. 2003). 

Each galaxy in this sample has been observed with the 
Spitzer IRAC and MIPS bands, in addition to observa¬ 
tions using the blue filter of the IRS peak-up facility. 
These blue filter peak-ups have spectral coverage from 
13.3 — 18.7 pm and give an additional photometric point 
at 16 pm, between the IRAC 8 pm and MIPS 24 pm 
bands. Aperture photometry was performed in 7" and 
12" radius apertures for the 3.6 — 8 ^m IRAC and 24 pm 
MIPS, respectively, and then aperture-corrected to 12.2" 
and > 35", respectively. For the MIPS 70 and 160 pm 
bands, nearly all the galaxies can be treated as point 
sources, and aperture corrections were taken from the 
MIPS handbook. A full description of these aperture 


corrections is described in Johnson et al. (2007). For a 
more detailed description of the SSGSS dataset, we refer 
the reader to O’Dowd et al. (2011). 

The IRS spectroscopy for the SSGSS sample is 
obtained through the NASA/IPAC Infrared Science 
Archive (IRSA) website® and is from the work of O’Dowd 
et al. (2011). This study utilizes the lower resolution 
Short-Low (SL) and Long-Low (LL) IRS modules, as 
these have been obtained for the entire SSGSS sam¬ 
ple. The SL module spans 5.2 — 14.5 ^m with resolving 
power R = 60 — 125 and has a slit width of 3.6 — 3.7". 
The LL module spans 14 — 38 ^m with resolving power 
R = 57 — 126 and has a slit width of 10.5 — 10.7". The 
spectra from the two modules were combined by weighted 
mean and a detailed description of the method can be 
found in O’Dowd et al. (2011). At z ~ 0.1 these galaxies 
are sufficiently distant such that the IRS slit encompasses 
a significant fraction of each galaxy (r-band Petrosian di¬ 
ameters are ^10"), providing some of the best MIR SEDs 
for a continuous SFR(A) determination. 

In order to accurately determine a diagnostic for star 
formation across the MIR, it is necessary to only consider 
cases where the majority of light is being contributed 
from stars (i.e., star-forming galaxies; SFGs) and not 
from an active galactic nucleus (AGN). The galaxy type 
is traditionally determined according to their location on 
the Baldwin-Phillips-Terlevich (BPT) diagram (Baldwin 
et al. 1981; Kewley et al. 2001; Kauffmann et al. 2003). 
By adopting the DR7 values of emission line measure¬ 
ments, we find that 64 of the 101 SSGSS galaxies are 
classified as SFGs. We note that the initial classification 
of SSGSS galaxies by O’Dowd et al. (2011) utilized SDSS 
DR4 measurements, which results in a few BPT designa¬ 
tions to differ between these works. We further exclude 
6 of the SSGSS galaxies classified as SFGs from our anal¬ 
ysis for the following reasons: SSGSS 18 appears to be a 
merger, SSGSS 19, 22, and 96 have significant breaks in 
their IRS spectra due to low signal-to-noise (S/N), and 
SSGSS 35 and 51 suffer from problems with IRS confu¬ 
sion. This leaves 58 galaxies to be used in our calibration 
of a monochromatic MIR SFR indicator. Table 1 shows 
the SSGSS IDs of the galaxies used for this study along 
with some of their properties (additional parameters in 
the table are introduced in later sections). 

Our sample of 58 galaxies span a redshift range of 
0.03 < z < 0.22 with a median redshift of 0.075. The 
range of infrared luminosity is 9.53 < log(AiR/A 0 ) < 
11.37, with a median of 10.55. All measurements of 
Air for these galaxies are taken from the original SS¬ 
GSS dataset (presented in Treyer et al. 2010). The selec¬ 
tion criteria for the SSGSS sample was based on 5.8 pm 
surface brightness and 24 pm flux density, and this re¬ 
stricts our sample of 58 galaxies to relatively high stellar 
masses (1.6 x 10® < M/Mq < 1.7 x 10^^) and metallic- 
ities (8.7 < 12 -|- log(0/H) < 9.2). These stellar mass 
and metallicity estimates are updated from the SSGSS 
dataset values (based on DR4) to the MPA-JHU DR7 
estimates. 

2.2. WISE Data 

The WISE All-Sky Survey provides photometry at 3.4, 
4.6, 12, and 22 pm (Wright et al. 2010) which comple- 


http: //www. stsci. edu/jwst/instruments/miri/instrumentdesigii/f liters/ 
http://www.mpa-garching.mpg.de/SDSS/DR7/ ® http://irsa.ipac.caltech.edu/data/SPITZER/SSGSS/ 


Continuous Mid-Infrared Star Formation Rate Indicators 


3 


TABLE 1 

Summary of Galaxy Properties and IRS Correction Terms 


SSGSS 

ID 

R.A. 

(J2000) 

Decl. 

(J2000) 

z 

log(LiR) 

(Lq) 

(SFR) 

(Mq/jt) 

*^phot 

k 

(16+fc) 

(8+fc) 

1 

160.34398 

58.89201 

0.066 

10.55 

4.22 

0.963 

14.10 

1.362 

2 

159.86748 

58.79165 

0.045 

9.83 

0.65 

1.062 

23.64 

1.253 

3 

162.41000 

59.58426 

0.117 

10.59 

3.41 

1.091 

54.09 

1.129 

4 

162.54131 

59.50806 

0.066 

10.22 

1.29 

0.871 

-4.466 

3.264 

5 

162.36443 

59.54812 

0.217 

11.37 

19.55 

1.088 

24.48 

1.246 

6 

162.52991 

59.54828 

0.115 

10.99 

7.41 

0.851 

1.310 

1.859 

8 

161.48123 

59.15443 

0.044 

9.97 

0.79 

0.843 

-5.224 

3.882 

14 

161.92709 

56.31395 

0.153 

11.06 

9.60 

0.918 

10.69 

1.428 

16 

162.04231 

56.38041 

0.072 

10.42 

2.15 

0.943 

15.73 

1.337 

17 

161.76901 

56.34029 

0.047 

10.83 

5.58 

0.922 

21.94 

1.267 

24 

163.53931 

56.82104 

0.046 

10.59 

3.43 

1.077 

0.240 

1.971 

25 

158.22482 

58.10917 

0.073 

10.30 

1.89 

1.061 

2.363 

1.772 

27 

159.34668 

57.52069 

0.072 

11.01 

7.70 

0.950 

-0.933 

2.132 

30 

159.73558 

57.26361 

0.046 

10.11 

1.06 

1.040 

63.01 

1.113 

32 

161.48724 

57.45520 

0.117 

10.82 

6.30 

1.121 

278.9 

1.028 

34 

160.30701 

57.08246 

0.046 

9.96 

0.79 

0.987 

46.47 

1.147 

36 

159.98523 

57.40522 

0.072 

10.39 

2.01 

0.996 

-1.629 

2.256 

38 

160.20963 

57.39475 

0.118 

10.92 

6.43 

0.905 

1.989 

1.801 

39 

159.38356 

57.38491 

0.074 

10.14 

1.28 

0.684 

1.615 

1.832 

41 

158.99098 

57.41671 

0.102 

10.34 

1.68 

0.818 

2.159 

1.787 

42 

158.97563 

58.31007 

0.155 

11.03 

9.32 

1.170 

14.84 

1.350 

46 

159.02698 

57.78402 

0.044 

10.02 

0.99 

0.932 

3.183 

1.715 

47 

159.22287 

57.91185 

0.102 

10.68 

4.34 

1.331 

61.86 

1.115 

48 

159.98817 

58.65948 

0.200 

11.24 

15.53 

0.869 

11.59 

1.408 

49 

159.51942 

58.04882 

0.091 

10.49 

3.13 

1.107 

0.495 

1.942 

52 

160.54201 

58.66098 

0.031 

9.53 

0.29 

1.129 

-2.795 

2.537 

54 

160.41264 

58.58743 

0.115 

11.20 

11.53 

0.901 

6.117 

1.567 

55 

160.29353 

58.25641 

0.121 

10.54 

2.95 

0.865 

3.155 

1.717 

56 

160.41617 

58.31722 

0.072 

10.01 

0.85 

0.886 

4.645 

1.633 

57 

160.12233 

58.16783 

0.073 

9.92 

0.75 

0.689 

15.28 

1.344 

59 

159.89861 

57.98557 

0.075 

10.36 

2.00 

0.928 

15.09 

1.346 

60 

160.51027 

57.89706 

0.116 

10.48 

2.89 

0.910 

6.530 

1.551 

62 

160.91280 

58.04736 

0.133 

11.08 

9.85 

0.978 

0.688 

1.921 

64 

161.00317 

58.76030 

0.073 

10.88 

5.44 

0.913 

3.276 

1.709 

65 

161.37666 

58.20886 

0.118 

11.18 

11.85 

0.965 

12.55 

1.389 

66 

161.25533 

57.77575 

0.113 

10.86 

6.75 

0.937 

13.86 

1.366 

67 

161.18829 

58.45495 

0.031 

10.09 

1.21 

1.581 

8.790 

1.476 

68 

163.63458 

57.15902 

0.068 

10.54 

3.37 

0.975 

25.56 

1.238 

70 

163.17673 

57.32074 

0.090 

10.49 

2.81 

1.037 

255.8 

1.030 

71 

163.21991 

57.13160 

0.163 

10.98 

8.32 

1.271 

285.0 

1.027 

72 

163.25565 

57.09528 

0.080 

10.79 

5.01 

0.951 

14.00 

1.364 

74 

161.95050 

57.57723 

0.118 

10.80 

5.07 

0.878 

-2.430 

2.436 

76 

162.02142 

57.81512 

0.074 

10.55 

2.62 

0.938 

6.128 

1.566 

77 

162.10524 

57.66665 

0.044 

9.69 

0.62 

1.130 

66.79 

1.107 

78 

162.12204 

57.89890 

0.074 

10.56 

3.17 

1.027 

-4.810 

3.508 

79 

161.25693 

57.66116 

0.045 

9.82 

0.81 

1.023 

-0.238 

2.031 

80 

162.07401 

57.40280 

0.075 

10.35 

2.14 

0.992 

-1.949 

2.322 

81 

162.04674 

57.40856 

0.075 

10.24 

1.69 

0.909 

1.709 

1.824 

82 

161.03609 

57.86136 

0.121 

10.77 

4.98 

0.854 

2.444 

1.766 

83 

160.77402 

58.69774 

0.119 

10.92 

6.56 

0.919 

-3.474 

2.768 

88 

161.38522 

58.50156 

0.116 

10.51 

2.54 

1.106 

42.43 

1.159 

90 

162.64168 

59.37266 

0.153 

11.02 

8.66 

0.796 

8.033 

1.499 

91 

162.53705 

58.92866 

0.117 

10.78 

5.31 

0.782 

-3.869 

2.937 

92 

162.65512 

59.09582 

0.032 

10.18 

1.42 

0.966 

6.675 

1.545 

94 

161.80573 

58.17759 

0.061 

10.13 

0.90 

0.919 

9.281 

1.463 

95 

163.71245 

58.39082 

0.115 

10.82 

5.74 

0.956 

22.68 

1.261 

98 

164.14571 

58.79676 

0.050 

10.59 

3.70 

0.982 

5.797 

1.580 

99 

164.33247 

57.95170 

0.077 

10.82 

5.58 

0.903 

-1.944 

2.321 


Notes. Columns list the (1) galaxy ID number, (2) redshift, (3) integrated infrared luminosity from 8 — 1000 /rm, (4) average SFR from 
the diagnostics in Table 2, (5) offset between IRS spectra and global photometry above 16 pm (6) correction parameter for wavelength- 
dependent aperture loss of IRS spectrum below 16 pm, (7) correction factor of the spectrum at 8 pm due to wavelength-dependent aperture 
loss. 
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ments the wealth of IR data available for the SSGSS sam¬ 
ple. Most importantly for this study, the WISE 12 /rm 
band provides a crucial photometric point that bridges 
the gap between the Spitzer 8 ^m and 24 ^m bands, 
a section of the MIR SED that experiences a transition 
from polycyclic aromatic hydrocarbon (PAH) emission 
features and dust continuum. The WISE photometry for 
the SSGSS sample is obtained through the NASA/IPAC 
Infrared Science Archive (IRSA) website^. Using an ap¬ 
proach similar to that of Johnson et al. (2007) for the 
Spitzer bands, we utilize the 13.75" radius aperture mea¬ 
surements for 12 /rm and then apply an aperture cor¬ 
rection of 1.20. This correction term was found using 
sources in our sample with no obvious contamination 
from neighbors and measuring the flux density out to 
24.75" to determine their total flux density. The WISE 
photometry at 22 ^m is less accurate than the Spitzer 
24 /im, owing to it having two orders of magnitude lower 
sensitivity (Dole et al. 2004; Wright et al. 2010), and is 
not used for our analysis. In addition, the 22 /rm ob¬ 
servations suffer from a effective wavelength error, which 
systematically brightens the photometry of star forming 
galaxies (Wright et al. 2010; Brown et al. 2014). 

3. ANALYSIS 


in this longer wavelength module, coupled with the larger 
slit width of 10.6" allows for the standard SPIGE algo¬ 
rithm to accurately recover the total flux density, since 
objects are close to point sources (see O’Dowd et al. 
2011). This would suggest that any offsets between the 
photometry and spectroscopy beyond 16 /iin should be 
uniform across the module (i.e., a global loss in flux den¬ 
sity). Therefore, each spectrum is first fit to match the 
16 pm and 24 pm photometric points using a constant 
offset, Cphot, found using chi-squared minimization. 


y^..-('S'nhot.?; *S'tRS.?,)/U'(* 5'nhot.7;) 




E 


*5’phot,z Cphot*5’lRS,^ 
^(•^photy) 


( 1 ) 

( 2 ) 


where S'photy is the Spitzer photometric flux density of 
band i, cr(S'phot,i) is the uncertainty of the Spitzer pho¬ 
tometric flux density, and iSiRs.i is the effective IRS pho¬ 
tometric flux density found using the transmission curve 
for each band, T'i(A), 


'S’lRsy 


/ ^iRs(A)r,(A)dA 
f T,(A) dX 


( 3 ) 


3.1. Anchoring the IRS Spectra to Global Photometry 

For this study, we focus on utilizing the Spitzer 5.8 — 
24 ^m and WISE 12 pm photometry to anchor the 
Spitzer IRS 5 — 40 ^m spectroscopy available for the 
entire sample. The reasoning for this approach is to 
have spectroscopy that is representative of the global 
flux density of each galaxy, which is required to create a 
calibrated continuous, monochromatic SFR(A) indicator. 
Offsets between the IRS spectrum and the photometry 
can occur from differences in data reduction methods, or 
from the width of the IRS slit being smaller than the size 
of the galaxy. The former effect results in a uniform off¬ 
set across the entire spectra and can be corrected with 
a normalization factor. The behavior of the latter ef¬ 
fect will be dependent on whether the galaxy observed 
is an unresolved point-like source. The default Spitzer 
IRS custom extraction (SPIGE) does include a correc¬ 
tion for light lost from the slit due to the changing an¬ 
gular resolution as a function of wavelength but assumes 
the object to be a point source. To correct for both of 
these effects, we utilize photometry from the Spitzer 8, 
16, and 24 pm and WISE 12 pm bands as a reference. 
The end-of-channel transmission drop of the SL mod¬ 
ule below ^^5.8 pm, combined with a typical redshift of 
z ^ 0.1, makes the 5.8 pm band region unreliable for 
use in most cases, and so it is not used as an anchor. 
However, we do make use of the 5.8 pm band to inspect 
our photometric matching in the lowest redshift galaxies 
(see below). Here we outline our approach to correct for 
offset effects so that these spectra are well representative 
of global photometric measurements. 

O’Dowd et al. (2011) found that IRS measurements of 
SSGSS galaxies from the LL module did not show evi¬ 
dence for significant aperture loss when compared to the 
Spitzer 16 and 24 pm photometry. This is attributed to 
the fact that the lower resolution (larger PSF) of sources 

^ http://irsa.ipac.caltech.edu/Missions/wise.html 


This method ignores the method of calibration that was 
utilized for each specific bandpass (i.e., a conversion of 
number of electrons measured by the detector into a flux 
density in Jy requires knowing the shape of the incoming 
flux of an object, which varies as a function of wave¬ 
length). However, discrepancies between the adopted 
method and correcting for calibration effects amount to 
^1%, and is not important for this study. The offset 
required to match photometric values is typically small, 
with values of Cphot being between 0.7 — 1.6. 

In contrast to the LL module, O’Dowd et al. (2011) 
found that IRS measurements short-ward of 16 pm from 
the SL module did show evidence for aperture loss when 
compared to the Spitzer 8 pm photometry. In this case, 
the increasing resolution of the SL module at shorter 
wavelengths results in many of the galaxies in this sample 
being resolved in this module. Also taking into account 
that the SL slit is 3.6", which is smaller than the average 
extent of ^10" (r-band Petrosian diameter) for SSGSS 
galaxies, implies that flux density loss in this wavelength 
region is more pronounced for more extended objects. 
For this reason, an additional correction term must be 
introduced below 16 pm which has a 1/A dependency 
to reflect the additional losses as resolution increases at 
shorter wavelengths (i.e., the PSF is decreasing at shorter 
wavelengths, resulting in less correction of light outside 
the slit). The correction terms adopted are summarized 
in the following equations, 

5irs,o.(A) = I X {W) : A < 16 /.m 

i >S'iRs(A)/cphot :\>lQpm 

(4) 

where fc is a constant found by performing a Levenberg- 
Marquardt least-squares fit of this function, using the 
IDL code MPFITFUN (Markwardt 2009), such that the IRS 
spectrum matches the 8 and 12 ^m photometric flux den¬ 
sity. Smaller values of k correspond to larger correction 
factors in the spectrum. Examples of normalizing the 
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IRS spectroscopy to the photometry are shown in Fig¬ 
ure I . A list of the normalization parameters is shown in 
Table 1. 

As a consistency check on the 1 /A dependency, a more 
accurate check is made on the few cases where spectra 
have high S/N and low redshifts, such that the 5.8 fim 
band region of the spectrum is reliable to use for convolu¬ 
tion. In nearly all these cases the correction using a 1/A 
dependency matches the observed photometric point, as 
demonstrated by SSGSS 46 in Figure 1. In addition, the 
agreement of the Spitzer 8 /rm and WISE 12 /im data 
with this approach suggests these normalized spectra are 
well representative of photometric values. 

3.2. Extending out to 70 fim 

The IRS spectrum extends out to 40 ^m, however, the 
combination of end-of-channel transmission drop around 
Aobs 37 /im and redshift effects makes these spectra 
unreliable for Arest ^ 37/(1 -I- z) /rm. For the mean red- 
shift of this sample (2 ~ 0 . 1 ), this corresponds to roughly 
Arest 34 /im. In order to utilize the wavelength region 
between 34 and 70 /im, which contains poor quality or 
no spectral data, we interpolate with dust models. This 
interpolation is expected to be robust since there are no 
sharp emission features in this wavelength range. The 
shape of the emission in this region is dependent on the 
temperature distribution of of the dust, the grain size dis¬ 
tribution, as well as the relative importance of stochastic 
versus thermal equilibrium heating (the former gives an 
almost-constant continuum and the latter is responsible 
to the Wien-like rise of the spectrum). 

To extend the wavelength region of our study out to 
70 /rm, we fit the dust models of Draine & Li (2007), 
combined with an additional stellar continuum compo¬ 
nent, to our IR photometry and IRS spectroscopy. For 
these models, the emission spectrum is given by Draine 
et al. (2007) as, 

5.,„,odel = L!*B.(r,) + [(1 - 

T M : Cniin; fAiaxj O^)? (5) 

where 11* is the solid angle subtended by stars, T* is the 
effective temperature of the stellar contribution, Mdust is 
the total dust mass, Dium is the distance to the galaxy, 
Pi, is the specific power per unit dust mass, C/min {Um&x) 
is the minimum (maximum) interstellar radiation, 7 is 
the fraction of the dust mass exposed to radiation with 
intensity U > Umin, Jm corresponds to the dust model 
(i.e., the PAH abundance relative to dust, (/pah; shown 
in table 3 of Draine & Li 2007), and a is the power- 
law factor for the starlight intensity. In summary, this 
emission spectrum is a linear combination of three com¬ 
ponents: ( 1 ) a stellar continuum with effective temper¬ 
ature T* which dominates at A < 5 /rm; (2) a diffuse 
ISM component with an intensity factor U = C/min; and 
(3) a component arising from photo-dissociation regions 
(PDRs). Typically, component (2) comprises a much 
larger amount of the total dust mass and, as such, is 
dominant over component (3) in the emission spectrum 
(Draine et al. 2007, 2014). 

To fit this model we follow the approach outlined in 
Draine et al. (2007), which found that the SEDs of 
galaxies in the Spitzer Infrared Nearby Galaxies Survey 


(SINGS) were well reproduced with fixed values of a = 2, 
Umax = 10®, and T* = 5000 K. Holding these parameters 
fixed, (jpAH, C/min, 7, Afdust, and H* are varied to find 
the dust model that comes closest to reproducing the 
photometry and spectroscopy. For this work, a grid of 7 
values is constructed for all gpAH (MW, LMC, SMG) and 
C/min values. The value of Mdust for each grid point is 
determined by minimizing the parameter in a similar 
manner to eq. ( 1 ) and ( 2 ), as Mdust represents a constant 
offset value. The goodness-of-fit for each case is assessed 
using the parameter, 

2 X ^ ‘^'obs.t *5'model,i /^\ 

X ^ L ^2 ^^2 - > (6) 

2 obs,i model,i 

where the sum is over observed bands and spectroscopic 
channels, Smodei.t is the model spectrum (for band com¬ 
parison, the model spectrum is convolved with the re¬ 
sponse function of that band), (Tobs,i is the observa¬ 
tional uncertainty in the observed flux density S'obs.i, 
and (Tmodei.i = O.lS'modei.i as adopted by Draine et al. 
(2007). The observed flux densities used in determining 
the best fit is comprised of the IRS spectroscopy in ad¬ 
dition to the IRAG 3.6, 4.5, and 5.8 pm and MIPS 70, 
and 160 pm photometry. The model which minimizes 
the value of x^ is adopted for use in representing the re¬ 
gion Arest ^ 37/(1 -I- 2 ) pm. Examples of the best fitting 
model for SSGSS galaxies are shown in Figure 2. The 
bands within the IRS region are shown only for compar¬ 
ison and are not used directly for the fit. 

Our choice to adopt the model which minimizes the 
value of is not necessarily the most accurate repre¬ 
sentation of the spectra, as the degeneracy of the model 
parameters can allow for multiple fits to have similar 
values while having different FIR SED shapes. However, 
we do not consider this to be of great significance for 
this study for two reasons. First, the flux density vari¬ 
ation due to changes in the SED shape for cases with 
ix^ ~ Amin) < where Xmin minimum value of the 

reduced is typically less than 25% over the 30 — 70 pm 
region, which is lower than the scatter among individual 
galaxy templates. As we will be utilizing an average of 
our galaxy templates for our diagnostic, the uncertainty 
from model SED variations will not be the dominant 
source of uncertainty. Second, the parameters of these 
fits are not used to determine the properties of these 
galaxies, which are more sensitive to these degeneracy 
effects than the total flux density. 

All cases are best fit by Milky Way dust models with 
(/PAH > 2.50%. There is a systematic trend among most 
fits to underestimate the flux density around the 8 pm 
PAH feature and overestimate the 10 — 20 pm region. 
This is most likely due to the limitations of fitting only 
three components to the data. However, since the main 
focus of these fits is to provide a description of the region 
from ~34 — 70 pm, these deviations are not considered to 
be significant, as they should have little effect on match¬ 
ing the shape of the emission beyond ~34 pm. As will 
be discussed in § 6 . 2 , these fits are consistent with other 
2 = 0 SFG templates found in the literature that make 
use of Draine & Li (2007) models, suggesting that these 
deviations could be a common problem. Investigation 
into the cause of these discrepancies warrants additional 
study, as it will improve our understanding of dust prop- 
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SSGSS1 z=0.066 SSGSS14 z=0.153 




SSGSS 46 z=0.044 


SSGSS 59 z=0.075 




^obs ^obs 

Fig. 1. — IRS spectroscopy (black line), Spitzer photometry (green squares), and WISE photometry (cyan circle) for some SSGSS galaxies. 
The IRS spectrum is normalized to the Spitzer 8 pm, 16 pm, and 24 pm and WISE 12 pm photometric flux densities according to the 
method described in § 3.1 (red line). The effective IRS photometry, found using the transmission curve for each bandpass filter, are 
shown as triangles. The normalized transmission curves for the Spitzer and WISE bands in this region are shown as green and cyan lines, 
respectively. 


erties in galaxies. We do not make use of these fits to de¬ 
termine Lir, values, instead using the values provided in 
the SSGSS catalog, which have been extensively checked 
(Treyer et al. 2010) across a variety of diagnostic meth¬ 
ods. 

3.3. Determining Rest-Frame Luminosities 

As these galaxies span a redshift range of 0.03 < z < 
0.22, the IRS spectrum and photometry of each galaxy 
span slightly different regions in rest-frame wavelength. 
This offset causes observed photometric values to vary by 
up to 30% from the rest-frame values. This would affect 
our determination of SFR if not accounted for and intro¬ 
duce additional scatter. Since previous MIR calibrations 
have been performed for local samples of galaxies {z ^ 0) 
to accuracies around 30% (Rieke et al. 2009; Kennicutt 
et al. 2009; Calzetti et al. 2010; Hao et al. 2011), this is 
a non-negligible effect. 

To correct for redshift effects, photometric values for 
each band are determined by convolving the spectrum 
at the rest-frame filter postions for each band accord¬ 


ing to eq (3), only now using the corrected spectrum, 
5'iRS,corr(A), instead of the original IRS spectrum. This 
is performed for the Spitzer 8, 24, and 70 /rm and WISE 
12 and 22 /rm bands. These corrected flux densities are 
used to calculate the rest-frame luminosity (erg s“^) of 
each band, 

Epest — (^7/p/)rest — (i^>5'lRS,corr)obs , (7) 

where Vrest and Vobs are the effective rest-frame and 
observer-frame frequency of each band, respectively, and 
Dium is the luminosity distance for the galaxy, calcu¬ 
lated from its redshift. These rest-frame luminosities are 
used to determine SFRs for each of the galaxies in our 
sample. In a similar manner, each IRS spectrum is ex¬ 
pressed as a wavelength dependent rest-frame luminos¬ 
ity, L(A)rest, using the continuous spectrum, S'iRS,corr(A). 
This is used later for calibrating our wavelength contin¬ 
uous SFR-luminosity conversion factors, C(A). 

To correct the Spitzer 3.6 and 4.5 /im bands, which 
lie outside of the IRS spectral coverage, a correction is 
applied assuming these bands encompass the Rayleigh- 


























































Continuous Mid-Infrared Star Formation Rate Indicators 


7 


SSGSS1 SSGSS14 




2 5 10 20 50 100 200 

A,est (“m) 


2 5 10 20 50 100 200 

Ares, (“m) 


SSGSS 46 


SSGSS 59 




2 5 10 20 50 100 200 

Ares, (“m) 


2 5 10 20 50 100 200 

Ares, (“m) 


Fig. 2.— A fit to the spectra of some SSGSS galaxies using the dust models of Draine Sz Li (2007), shown as the solid red line. The IRAC 
3.6, 4.5, and 5.8 fim along with the MIPS 70 and 160 /rm photometry is used to fit the dust continuum in the absence of the IRS spectrum. 
The bands within the IRS region are shown for comparison and not used directly for the fit. The regions of the IRS spectrum associated 
with transmission drops in the instrument are not shown for clarity. We attribute the relatively poor match in the 6 — 20 fim region to 
using a simple three component model. However, the purpose of these fits is only to determine the shape of emission in the 34 — 70 fim 
region, for which the data is found to be in good agreement with the models. 


Jeans tail of the stellar continuum emission, 

= ^obs (= 5obs X (1 + . (8) 

\ Xobs / 

The luminosity is then found following eq (7) using 
The rest-frame Spitzer 3.6 and 4.5 /im luminosities will 
only be used to examine the origins of scatter within 
our conversion factors in § 6.1, and has no influence on 
our estimates for the MIR conversion factors. The fit¬ 
ting results in § 3.2 suggest that this simple approach is 
reasonable for our sample of SFGs. 

3.4. Reference Monochromatic SFR Indicators 

In order to perform any calibration of luminosity as 
a SFR indicator, it is necessary to rely on previous, 
well-calibrated SFR indicators. In this work, we uti¬ 
lize the calibrations of Kennicutt et al. (2009), Rieke et 
al. (2009), Calzetti et al. (2010), and Hao et al. (2011), 
which incorporate the full suite of data available for this 


sample. This list is shown in Table 2. The reference 
SFR, (SFR), for each galaxy is taken to be the average 
of the SFRs from these calibrations. By utilizing the av¬ 
erage of a large number of diagnostics, we limit the risk 
of potential biases that any single diagnostic can be sub¬ 
ject to. Several of these reference diagnostics make use 
of the total infrared luminosity, Ttir, which refers to the 
integrated luminosity over the region from 3 to 1100 p-m. 
All measurements of Ttir for these galaxies have been 
obtained from the original SSGSS dataset (Treyer et al. 
2010). For the SSGSS sample, Treyer et al. (2010) find 
that Ltir is larger than Lir by ^0.04 dex. 

To utilize SDSS measurements of Ha for a SFR esti¬ 
mation, it is necessary to apply an aperture correction. 
The diameter of the SDSS spectroscopic fiber spans 3", 
which is a factor of '^3 smaller than the typical size of 
the SSGSS galaxies and results in only a fraction of the 
light being measured. We correct for this aperture effect 
using the prescription from Hopkins et al. (2003), which 
uses the difference between the r-band Petrosian mag- 













































Battisti et al. 


TABLE 2 

Reference Star Formation Rate Calibrations 


Band(s) 

Lx Range 
(erg s-l) 

SFR 

logCj, 

Reference 

FUV-I-TIR 


[F(FUV)obs -b 0.46L(TIR)]/C,, 

43.35 

Hao et al. (2011) 

FUV+24 iim 


[i(FUV)obs-b3.89L(24)]/C, 

43.35 

Hao et al. (2011) 

NUV-I-TIR 


[L(NUV)obs -b 0.27L(TIR)]/C,, 

43.17 

Hao et al. (2011) 

NUV-l-24 /LHi 


[t'(NUV)obs + 2.26L(24)]/C, 

43.17 

Hao et al. (2011) 

HqH- 8 /im 


HlHalobs -b 0.011 LCSll/C,, 

41.27 

Kennicutt et al. (2009); Hao et al. (2011) 

Ha-l-24 iim 


[i(Ha)obs +0.020 L(24)]/C,i 

41.27 

Kennicutt et al. (2009); Hao et al. (2011) 

Ha-l-24 iim 

L(24) < 4 X 

[i(Ha)obs-b 0.020 L(24)]/C,, 

41.26 

Calzetti et al. (2010) 


4 X 10"^^ < L(24) < 5 X 

[i(HQ)obs +0.031 L(24)]/C,, 

41.26 

Calzetti et al. (2010) 


L(24) > 5 X lO"*® 

L(24) X [2.03 X 10-‘‘‘L(24)]0 “'‘8/c„ 

42.77 

Calzetti et al. (2010) 

Ho-I-TIR 


[L(H«)obs + 0.0024 L(TIR)]/C:, 

41.27 

Kennicutt et al. (2009); Hao et al. (2011) 

24 urn 

2.3 X 10'‘2 < L(24) < 5 X 

L{2A)ICx 

42.69 

Rieke et al. (2009) 


L(24) > 5 X lO"*® 

L(24) X (2.03 X 10-'‘'‘L(24))° °48/c'„, 

42.69 

Rieke et al. (2009) 

70 pm 

L(70) > 1.4 X 

L(70)ICx 

43.23 

Calzetti et al. (2010) 


Notes. Columns list the (1) bands used in the calibration, (2) luminosity range over which the calibration can be used; empty fields denote 
an unspecified range, (3) SFR conversion formula, (4) conversion constant, and (5) reference for calibration. 


nitude and the r-band fiber magnitude (see also Treyer 
et al. 2010). The values of these corrections range from 
1.9 — 8.3 for our sample, with the exception of SSGSS 67 
with a correction of 21.4 due to its much larger size. 

As a consistency check, the SFRs inferred from each 
indicator for our galaxies is compared in Figure 3. It is 
seen from the distribution that the majority of these val¬ 
ues agree within the ^^30% uncertainty associated with 
individual calibrations (Rieke et al. 2009; Kennicutt et 
al. 2009; Calzetti et al. 2010; Hao et al. 2011). A formal 
fit of this distribution to a Gaussian profile gives values 
oi fj, = —0.02 and tr = 0.17. Our choice in using the 
average of the diagnostics, {SFR), for each galaxy in¬ 
stead of the median value appears to cause no significant 
differences, with typical offsets of only a few percent be¬ 
tween the two, which are symmetric. A formal Gaussian 
fit of fractional difference between the mean and median 
gives R = —0.01 and a = 0.04. We illustrate the rel¬ 
ative offsets for the individual calibrations in Figure 4. 
We note that the MPA/JHU group provides independent 
estimates for SFRs based on the technique discussed in 
Brinchmann et al. (2004) using extrapolated Ha mea¬ 
surements. However, we find that the spread in values 
for these estimates relative to (SFR) are significantly 
larger than our other diagnostics {la = 0.48), which we 
attribute to their larger SFR uncertainties (^50%), and 
as such were excluded from our analysis (including Fig¬ 
ures 3 and 4). 

3.5. Lir as a SFR Indicator 

A commonly utilized method to determine SFRs for 
galaxies relies on measuring the integrated luminosity 
over most of the IR wavelength range, Air(8 —1000 ^m). 
However, physically understanding the conversion of Lir 
to a SFR is non-trivial and sensitive to many assump¬ 
tions, such as the timescale of star formation, r, the star 
formation history (SFH), the metallicity, and the initial 
mass function (IMF; see Murphy et al. 2011; Calzetti 
2013). For example, a galaxy with a constant SFH, a 
fixed metallicity, and a fixed IMF will have the calibra¬ 
tion constant for Lir change by a factor of 1.75 between 
assuming r = 100 Myr and r = 10 Gyr (Calzetti 2013). 

We chose to avoid the use of SFRs based solely on 


Lir for reference because of the sensitivity to these as¬ 
sumptions. However, in order to compare the accuracy 
of our calibration on higher redshift samples (in § 6.4) 
for which Lir is the only technique available to esti¬ 
mate SFRs, we use a SFR-Lir conversion which repro¬ 
duces the values of (SFR) seen for the SSGSS sample. 
This occurs for a conversion factor of log[C(LiR)] = 
43.64 erg s“^/(M 0 yr“^). Utilizing Starburst99 (Lei- 
therer et al. 1999), with a constant SFH, solar metallicity, 
a Kroupa IMF over 0.1 — 100 Mq, and assuming all of 
the stellar light (UV-bvisible) is reradiated by dust, this 
corresponds to a timescale of r ~ 500 Myr (e.g., Calzetti 
2013). 

This adopted conversion factor differs slightly from 
other commonly adopted values. In the case of Murphy 
et al. (2011), log[C(LiR)] = 43.41 erg s“V(-^©yr”^), our 
calibration is larger by 70%. This large difference is due 
to two reasons: (1) Murphy et al. (2011) assume that 
only UV light is being reradiated by the dust and does 
not account for the optical light that would also be rera- 
dited (~40% of the discrepancy), and (2) they assume 
a 100 Myr constant star-forming population (^20% of 
the discrepancy). In the case of Kennicutt (1998) af¬ 
ter converting from a Salpeter (1955) IMF to a Kroupa 
(2001) IMF, log[C'(LiR)] = 43.53 erg s-^/{M qjt-^), our 
calibration is larger by 30%. Most of this difference is 
due to them assuming a 100 Myr constant star-forming 
population. 

4. A CALIBRATED CONTINUOUS, 
MONOCHROMATIC SFR(A) 

4.1. Composite IRS Spectrum 

The SFR of a galaxy, using a calibrated single-band 
luminosity, can be written as 

SFR(M 0 yr-i) = L^/C, , (9) 

where is the monochromatic luminosity, in units of 
erg s“^, and Cx is the conversion factor between SFR 
and luminosity for filter x (following convention of Ken¬ 
nicutt & Evans 2012). In this respect, the appropriate 
conversion factor at a given band is found by normaliz¬ 
ing the luminosity by the SFR determined independently 
from a reference calibration. 
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Fig. 3. — Left: Comparison of SFRs determined from the calibrations listed in Table 2 for the SFGs in the SSGSS. Each vertical strip 
of values shows the SFR values for each method on a single galaxy. The reference SFR, (SFR), for these galaxies is taken to be the 
average of the SFR values from these calibrations. Right: Histogram showing the distribution of SFR offsets relative to the reference value. 
This distribution is well fit by a Gaussian with p = —0.02 and cr = 0.17, suggesting that the majority of these values agree within the 
uncertainties associated with the individual calibrations. 


This same approach is taken to calibrate our continu¬ 
ous wavelength conversion factors, 

C'(A)(erg s-V(MQyr-i)) = r(A),est/(SFR), (10) 

where L(A)rest is the wavelength dependent IRS lumi¬ 
nosity and (SFR) is the reference SFR. To achieve our 
calibration of C'(A), the SFR-normalized IRS spectra are 
averaged together to create a composite spectrum for 
the group. As a result of shifting the spectra to the 
rest-frame, the wavelengths associated with each spectral 
channel no longer match exactly. Therefore, to perform 
this average, the channel wavelengths in the spectrum 
of the first galaxy in our group is taken to be the refer¬ 
ence grid. Next, the normalized luminosity values of the 
other galaxies are re-gridded to this (i.e., each channel is 
associated to the nearest neighboring reference channel). 
In using this approach, the smoothing of sharp features 
that result from direct interpolation is avoided. The un¬ 
certainties associated with this re-griding to determine 
a composite spectrum are small relative to the channel 
flux density uncertainly. Furthermore, these uncertain¬ 
ties are much smaller than the scatter between spectra, 
which drives the uncertainty of our template, and can be 
considered negligible for the purposes of this study. 

The result of an average for the entire sample of SFGs 


in the SSGSS sample is shown in Figure 5. The uncer¬ 
tainty of each channel in the composite spectrum is taken 
to be the standard deviation of the the group value for 
that channel. The sample standard deviation is the dom¬ 
inant source of uncertainty (typically between 20 — 30% 
of the normalized luminosity value) and is larger than 
the flux density uncertainties of individual IRS channels 
(typically ^2%) by roughly an order of magnitude. This 
template can be used to determine the appropriate con¬ 
version factor for any luminosity within our wavelength 
coverage. 

4.2. Filter Smoothed Composite Spectrum 

In practice, observations of a galaxy are made using 
specific bandpass filters that encompass a portion of their 
SED. Therefore, it is more practical to utilize a compos¬ 
ite spectrum that corresponds to photometric luminosi¬ 
ties observed by various bands as functions of redshift. 
To accomplish this, the normalized IRS spectrum of each 
SFG in our sample is convolved with the filter response of 
specific bands as functions of redshift (i.e., the effective 
wavelength blue-shifts and the bandpass narrows, both 
by a factor of (1+^), as one goes to higher redshifts). Per¬ 
forming this convolution is similar to smoothing by the 
bandpass filter, only with the filter width changing with 
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Fig. 4.— Histograms showing the distribution of the individual SFR estimates relative to the reference value, (SFR). The calibrations 
being considered are listed in Table 2. The median offset and Icr dispersion are shown in each panel. 


redshift. Throughout the rest of this paper, the term 
“smoothed” is used interchangeably to mean this convo¬ 
lution process. The redshift limit imposed for each band 
occurs at the shortest usable rest-frame wavelengths of 
the IRS spectrum for that band. The composite IRS 
template and the filter smoothed composites presented 
in this section are publicly available for download from 
the IRSA8. 

Each normalized IRS spectra is smoothed using the 
Spitzer^'^^ , WISE^^, and JWST/MIRI filters. The prop¬ 
erties of these filters is listed in Table 3. We note that the 
Herschel PACS 70 ^m is close enough to Spitzer 70 /rm 
that these can be interchanged for use with C' 57 o(A). We 
emphasize to the reader that care should be taken when 
considering the 22 /iin band, as it has been shown to suf¬ 
fer from an effective wavelength error (see Wright et al. 
2010; Brown et al. 2014). For the MIRI filters, we use 
the response functions of Glasse et al. (2015). Since these 
curves do not take into account the wavelength depen¬ 
dent quantum efficiency, the instrument transmission, or 
the responsivity of the detector, they should be updated 
once better curves become available. 

The composite spectra for each of these bands is cre- 

® http://irsa.ipac.caltech.edu/data/SPITZER/MIR_SFR 

® http;//irsa.ipac.caltech.edu/data/SPITZER/docs/irac/ 
http://irsa.ipac.caltech.edu/data/SPITZER/docs/mips/ 
http://www.astro.ucla.edu/-wright/WISE/passbands.html 


ated by averaging the smoothed spectra together, in the 
same manner as for the group average. The result of av¬ 
eraging the convolved spectra is shown in Figures 6 & 
7. The associated uncertainty with each channel (sam¬ 
ple standard deviation) is slightly lower than the native 
composite spectrum owing to the smoothing from the 
convolution and is typically between 15 — 20% of the nor¬ 
malized luminosity value (except for 70 pm case, which 
is still around 30%), making them comparable to accura¬ 
cies achieved in many previous calibrations. Previously 
determined MIR conversion factors (from z ^ 0 samples) 
are also shown and appear in good agreement. 

For the filter bands considered here, the smoothed 
IRS spectra show a very large increase in scatter below 
^6 pm, which is due to a combination of the end-of- 
channel uncertainties being very high and also from vari¬ 
ations in the old stellar populations of these galaxies. 
For these reasons, we only consider the regions for which 
the Icr uncertainty is less than 30% suitable for calibra¬ 
tion. In the case of the WISE 12 pm band, this region 
occurs below '^7 pm because of the significantly wider 
filter bandwidth. The ranges chosen for the calibration 
of each band is shown in Table 4. 

4.3. Fits to the Composite Speetra 

To simplify the application of our results as SFR in¬ 
dicators, each of the filter smoothed composite spectra 
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Fig. 5. — Normalized IRS luminosity, L(A)rest/{5'-FR), for all SFG galaxies (gray solid lines). The composite spectrum of this group is 
shown (thick black line) along with the standard deviation from this average (black dotted lines), which at most wavelengths is between 
20 — 30%. The fits to the dust continuum for each galaxy (gray dashed lines; described in § 3.2) along with the average (black dashed line) 
are also shown. The low dispersion among normalized spectra suggests that the 6 — 70 pm region can be utilized for SFR diagnostics. 


TABLE 3 
Filter Properties 


Instrument 

Band 

(^tm) 

FWHM 

(^tm) 

IRAC 

8 pm 

7.87 

2.8 

MIPS 

24 pm 

23.68 

5.3 

MIPS 

70 fj,m 

71.42 

19.0 

WISE 

12 pm 

12.08 

8.7 

WISE 

22 pm 

22.19“ 

3.5 

MIRI 

FIOOOW 

10.00 

2.0 

MIRI 

F1280W 

12.80 

2.4 

MIRI 

F1500W 

15.00 

3.0 

MIRI 

F1800W 

18.00 

3.0 

MIRI 

F2100W 

21.00 

5.0 

MIRI 

F2550W 

25.50 

4.0 


Notes. Columns list the (1) instrument, (2) passband name, (3) 
rest-frame effective wavelength, and (4) full width at half maxi¬ 
mum. “^The 22 pm observations suffer from a effective wavelength 
error (see Wright et al. 2010; Brown et al. 2014). 


is fit using a continuous function, fitx{X). We perform 


Levenberg-Marquardt least-squares fits of a polynomial 
(up to 1st order) and Drude profiles (up to 5), /r(A), 
to the smoothed composite spectra, using the IDL code 
MPFITFUN, 


2 5 

/zf,(A) = (II) 

i—1 r=l 


where x corresponds to the filter being considered, Pi 
are constants, and Ir{X) are Drude profiles. Drude pro¬ 
files, which are typically employed to characterize dust 
features, have the form 


T(A) 


(A/A, 


b 

W Ir _ 

-A,/A)2+72’ 


( 12 ) 


where A, is the central wavelength of the feature, 7 , is 
the fractional FWHM, and 6 , is the central intensity, 
which is required to be non-negative. We emphasize that 
because these are smoothed spectra, the parameters of 
these fits are not of physical significance and are simply 
being employed for ease of application. 
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For the Drude profiles, the central wavelengths, A^, 
are fixed to wavelengths that roughly correspond to the 
peaks in the smoothed spectrum, while jr and are 
left as free parameters. Therefore, there are up to 12 
free parameters in total, two from the polynomial and 
ten from the Drude profiles. The values of A^. for each 
smoothed composite fit and all the other fit parameters 
are listed in Table 4. The fits are shown for the individual 
bands in Figure 6 . The fitting functions are typically 
accurate to within ±5% (0.02 dex) of the true values and 
can be used in place of the templates [(^^(A) = fitx{X)]- 

4.4. Comparison to WISE SFR Calibrations 

The WISE All-Sky Survey (Wright et al. 2010) pro¬ 
vided photometry for over 563 million objects, and as 
such has great potential for future application of our cal¬ 
ibrations. Recently, calibrations of the WISE bands as 
SER indicators have emerged (Donoso et al. 2012; Shi 
et al. 2012; Jarrett et al. 2013; Lee et al. 2013; Cluver 
et al. 2014), some of which can be easily compared to 
our results. In particular, we focus on the results of Jar¬ 
rett et al. (2013) and Lee et al. (2013) as these have 
linear calibrations of the WISE 12 and 22 pm bands and 
can be directly compared to our values. The difference 
between the calibration values found in these works is 
rather large, corresponding to 0.34 dex (~120%) and 0.15 
dex (^40%), for the WISE 12 and 22 ^m band, respec¬ 
tively. These large discrepancies are likely the result of 
the different approaches of the two works. Jarrett et al. 
(2013) rely of the previous calibrations of Rieke et al. 
(2009) at 24 ^m, whereas Lee et al. (2013) attempt to 
determine SFRs from extinction-corrected Ha emission. 

The composite of our sample of SFG spectra smoothed 
by the WISE 12 and 22 pm filters is compared to these 
calibrations in Eigure 6 . We find that the results lie in- 
between the values found by Jarrett et al. (2013) and Lee 
et al. (2013). Since the WISE 22 pm band is so similar 
in shape and location to the Spitzer 24 ^m band, the 
calibrations of Zhu et al. (2008) and Rieke et al. (2009) 
are also presented and show close agreement to our work. 

5. APPLICATION TO HIGHER REDSHIET 
GALAXIES 

5.1. Demonstration 

Here we demonstrate how to apply our calibrations to a 
SEG with a known redshift. Let us consider using the ob¬ 
served 24 ^m flux density for the galaxy SSGSS 1 to esti¬ 
mate the SFR of this galaxy. At z = 0.066 for this partic¬ 
ular galaxy, the 24 ^m band has an effective wavelength 
of Aeff(24) = 22.51 pm and an observed flux density of 
5068(24) = 9.72 X 10“^® erg s“^ cm“^ Hz“^, which corre¬ 
sponds to a luminosity of log[Lobs(24)] = 43.11 erg s“^. 
Knowing the effective wavelength, we next want to use 
the composite 24 pm band smoothed spectrum to deter¬ 
mine the appropriate conversion factor, which is found 
to be log[C24(22.51 pm)] = 42.62 erg s“^/(M 0 yr“^) 
using the smoothed composite template or the fitting 
function (see Figure 6 ). Finally making use of eq. (9), 
we get that the SFR is simply the observed luminosity 
in this band divided by the conversion factor, SFR = 
1043.il/io42 .62 _ 3 09 M 0 yr“^. This value differs from 
the actual value of (SFR) = 4.22 M 0 yr“^ by about 30%. 

Next we can consider the slightly more distant case 
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Redshift Redshift 

0.31 0,12 0,0 2.95 1.96 1.37 0.97 0.69 0,48 0.320,18 0,0 




Redshift 


Redshift 




A^(S70) (urn) 


log[4,(W12) (um)] 


Redshift 


2.70 1.77 1.22 0.85 0.59 0.39 0.23 0.110.0 



Fig. 6. — The top of each panel shows the conversion factor for a Spitzer or WISE band using all SFG galaxies (solid red lines), along 
with their uncertainty (dotted red lines), which for most cases is between 15 — 20%. The solid black line is a fit to the smoothed spectrum, 
fitx(X). Local conversion factors from the literature are also shown for comparison (colored symbols). The region below -^6 pm is excluded 
due to significantly increased uncertainty in the composite spectrum (see § 4.2). The bottom of each panel shows the residuals between 
the conversion factor and a fit to the curve {\og[Cx{X)/fitx(X)]). 
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of SSGSS 14, at z = 0.153, and determine the SFR 
from its observed 24 ^m flux density of Sobs (24) = 
4.44 X 10“^® erg s“^ cm“^ Hz“^, which corresponds 
to a luminosity of log[Lobs(24)] = 43.55 erg s“^. The 
24 ^m band has an effective wavelength of Aeff(24) = 
20.53 ^m, which corresponds to log[C'24(20.53 /rm)] = 
42.57 erg s~^/(MQyr~^). Taking the ratio of these num¬ 
bers gives SFR = 9.63 MQyi~^. This value differs from 
the actual value of (SFR) = 9.60 M0yr“^ by < 1%. 

In the same manner, each calibration can be applied 
to any redshift that spans the Aeff range covered by IRS. 
These examples highlight the importance of using large 
sample sizes in the application of these diagnostics, as a 
single case can have SED variations relative to the mean 
of our sample, which can give rise to slight inaccuracies 
in SFR estimates. It is important to emphasize that 
the accuracy of such an application is dependent on the 
shape of the SED of SEGs as a function of redshift. The 
extent to which this condition holds is examined in detail 
in § 6.2. 

5.2. Limitations of this Sample 

It is important to acknowledge the potential differences 
of this sample with respect to high-z galaxies as well as 
the limitations for its use. As was mentioned, the selec¬ 
tion criteria for the SSGSS sample limits it to relatively 
high metallicities, which may not be a well representa¬ 
tive sample as one goes to high-z. In addition, if the 
dust content of high-z galaxies is different, it is possible 
that the amount of UV light reprocessed by dust could 
change. Eor example, if high-z galaxies had more dust, 
then our templates would overestimate the SFR, as it 
would be implicitly adding back in unobscured UV flux 
present in the SSGSS sample but that may not be there 
for the high-z galaxies. Variations in the typical dust 
temperature of galaxies with redshift would also pose a 
problem, as this would result in variations in their FIR 
SED. The relative importance of some of these effects will 
be tested when we compare our SED to those at higher 
redshift (§ 6.2). 

Another area for concern is in the range of Lm values 
spanned by the SSGSS sample. Our template is made 
utilizing galaxies over a range of 9.53 < log(LiR/L0) < 
11.37, which is lower than the range that is currently ac¬ 
cessible at high-z. However, the results of Elbaz et al. 
(2011) suggest that luminous infrared galaxies (LIRGs; 
Lir, > 1O^^L0) and ultra luminous infrared galaxies 
(ULIRGs; Lir > 1O^^L0) identified in the GOODS- 
Herschel sample at high-z have similar SEDs to nor¬ 
mal SEGs, in contrast to their starburst-like counter¬ 
parts found locally (e.g., Rieke et al. 2009). They find 
that the entire population of IR-bright galaxies has a 
distribution with a median of IRS = LiR/Lrest(8/rm) = 
4.9 [—2.2, -1-2.9], where the term in brackets is the Icr dis¬ 
persion. Elbaz et al. suggest that this population can be 
separated into two groups: a main-sequence (MS) of nor¬ 
mal SEGs for which IRS = 4±2 consisting of ~S0% of the 
population, and starburst (SB) galaxies which occupy 
the region with IRS > 8 and represent about ^20% of the 
population. Eor reference, the IRS value of our template 
is 4.8, which agrees with the median of the GOODS- 
Herschel sample. The uniformity of IRS values in MS 
galaxies, over the range 10® < Li^/Lq < 10^®, suggests 


that the SED of normal SEGs do not change drastically 
with luminosity. This also indicates that our limited 
range in Lir coverage for the SSGSS sample should not 
drastically affect its utility towards higher luminosity MS 
galaxies. 

Perhaps the biggest factor limiting the large scale ap¬ 
plication of this technique is in the ability to identify 
galaxy types at higher redshifts. The calibrations pre¬ 
sented in this work are applicable to normal SEGs, typ¬ 
ically referred to as being on the main-sequence of star 
formation, and not to cases undergoing starburst activ¬ 
ity (different SED) or with AGN (significant IR emission 
not associated with star formation). This topic should 
be thoroughly addressed before widespread applications 
of these calibrations can be made to specific surveys. 

There are a several techniques that have been demon¬ 
strated to isolate out AGN and starburst galaxies, how¬ 
ever, some of them rely on observations made outside 
the MIR. One of the most reliable techniques to identify 
AGN is thorough X-ray observations (e.g., Alexander et 
al. 2003), however these can miss obscured AGN and 
could be biased (Brandt & Hasinger 2005). In order to 
avoid obscuration effects, AGN selection techniques us¬ 
ing the MIR and FIR have also been developed. These 
include Spitzer+Herschel color-cuts (Kirkpatrick et al. 
2012), Spitzer/IKAG color-cuts (Lacy et al. 2004; Stern 
et al. 2005; Donley et al. 2012; Kirkpatrick et al. 2012), 
and WISE color-cuts (Stern et al. 2012; Mateos et al. 
2012; Assef et al. 2013). Emission line diagnostics, such 
as the BPT diagram (Kewley et al. 2013a) and the Mass- 
Excitation diagram (Juneau et al. 2014), are also effective 
techniques. It has been suggested that starburst galaxies 
can be identified as sources with IRS > 8 by Elbaz et al. 
( 2011 ). 

6 . DISGUSSION 

6.1. Origins of the Scatter in SFR(X) 

In addition to grouping all of the SEGs together, we 
also examine grouping our galaxies based on their lumi¬ 
nosity at rest-frame 3.6, 4.5, 8, 24, and 70 /rm, as well as 
their Lir, Lir surface brightness, and L(Hq!)/L( 24 ^m) 
ratios, in order to identify possible origins to the scatter 
within the SER calibration of the entire group. These 
bands are chosen because 3.6 ^m and 4.5 ^m correlate 
with the underlying stellar population (i.e., stellar mass; 
Meidt et al. 2012, 2014), and the other bands correlate 
strongly with star formation (Zhu et al. 2008; Rieke et al. 
2009; Calzetti et al. 2010; Hao et al. 2011). For each of 
these cases, the sample is divided into 6 bins with 9 — 10 
galaxies in each. 

Looking at each of the calibrations, weak trends are 
found suggesting larger conversion factors, at almost all 
MIR wavelengths, for galaxies with higher luminosities 
when arranged by any of the luminosities mentioned be¬ 
fore. A few examples are shown in Figure 8. These trends 
are very weak because the separation between the groups 
is comparable to the scatter within each of the groups, 
which is 10 — 25% (the largest scatter at lowest lumi¬ 
nosity galaxies), and similar to the uncertainty of the 
entire group average values. If real, these trends could 
suggest that (I) galaxies with a larger old stellar popula¬ 
tion require slightly larger conversion factors at all MIR 
wavelengths, as might be expected if light unassociated 
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Redshift 

0.67 0.25 0.0 



Redshift 
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-l^„(F1280W) («m) 
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Redshift 



5 6 7 8 9 10 20 30 


7,,(F2100W) iiim) 


Redshift 
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Fig. 7. — The top of each panel shows the conversion factor for select JWST/MIRI bands using all SFG galaxies (solid red lines), along 
with their uncertainty (dotted red lines), which for most cases is between 15 — 20%. The region below ~6 pm for each band is excluded 
due to significantly increased uncertainty in the composite spectrum (see § 4.2). The bottom of each panel shows the residuals between 
the conversion factor and a fit to the curve (log[C' 2 ;(A)// 2 tcc(A)]). 
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to star formation is contaminating the MIR; and/or (2) 
a slightly super-linear relationship exists between lumi¬ 
nosity and SFRs in the MIR. 

Our attempts to account these effects by introducing 
additional terms into the conversion factors do not ap¬ 
pear to signihcantly reduce the overall scatter of the con¬ 
version factors. Such additional terms would also make 
application of this method to higher redshift galaxies 
more difficult, as more information would be needed (e.g., 
determination of rest-frame luminosities). Given our lim¬ 
ited range in galaxy properties to determine the validity 
of any trends, we adopt the simplest approach and use 
the entire group average for our analysis. We refrain 
from using higher luminosity local galaxies, such as the 
(U)LIRGs in the GOALS sample (Armus et al. 2009), as 
an additional test of such claims because of the signifi¬ 
cant FIR SED evolution that occurs for Lir > lO^^L©, 
which is absent from galaxies of these luminosities at 
high-z (see § 5.2; Elbaz et al. 2011). In addition, a large 
fraction of these systems are likely to host AGN (U et al. 
2010 ), which would be excluded by our selection process. 
For reference to the reader, we illustrate the local SED 
evolution by comparing the GOALS photometry (U et 
al. 2010) and a few of the Rieke et al. (2009) templates 
to our own template, normalized by Lrest(8/rm), in Fig¬ 
ure 9. 

6.2. Variation in SFG SEDs with Redshift 

Many studies have sought to characterize the SED of 
different galaxy types (e.g., SFG, AGN) as functions of 
redshift. In this section, the templates of Elbaz et al. 
(2011); Kirkpatrick et al. (2012); Magdis et al. (2012) 
and Ciesla et al. (2014) are compared to our own to 
thoroughly examine the extent to which the SED of SFGs 
change with redshift. Similar to the approach outlined in 
this work, these studies use large surveys to construct IR 
templates for different populations of galaxies at different 
redshifts. In general, the templates created in these stud¬ 
ies suggest that the mean dust temperature of galaxies 
increases as one looks to higher redshifts. In addition to 
the change in dust temperature that is evident, Magdis et 
al. (2012) suggest that the value of IRS increases mildly 
from IRS ^ 4 to IRS ~ 6 at z > 2 for MS galaxies. 

If one considers the notion that both Lrest(8/rm) and 
Lir, are typically used for SFR indicators, such a change 
in IRS would suggest that there is a change in the SFR 
converstion factor of one (or both) of these luminosities 
with redshift and this is important to keep in mind when 
comparing the templates. Changes in Li.est(8/im) could 
result from variations in PAH abundances relative to the 
total dust content, which has been found to correlate 
with metallicity (Engelbracht et al. 2005, 200S; Marble 
et al. 2010), and also to the hardness of the radiation field 
(Madden et al. 2006; Gordon et al. 200S; Engelbracht et 
al. 2008). Given the sensitivity of Air to the contribu¬ 
tions from older stellar populations (Calzetti et al. 2010), 
it is also likely that the value of the conversion factor for 
SFR-Lir could also vary with redshift. 

The comparison between the templates from the lit¬ 
erature to our own is shown in Figures 10 and II. We 
have chosen to normalize the templates in two ways, both 
of which correlate with star formation. Normalizing by 
a close proxy for star formation is crucial to compare 
how viable our continuous calibrations are at higher red¬ 


shifts. The first method is to normalize by Lrest(8/rm), 
which is chosen over use of the 24 fj,m region because it is 
not available for the Kirkpatrick et al. (2012) templates. 
With this choice of normalization, it is also easier to di¬ 
rectly compare the shape of the MIR SEDs. The second 
method is to normalize by Air, as is traditionally done 
in many template comparisons. We reemphasize that the 
observed trend of IRS increasing from 4 to 6 implies that 
these choices of normalization for the templates will give 
different results. 

First, the templates of Kirkpatrick et al. (2012) are ex¬ 
amined as these provide the best sample for comparison 
because they are based on direct spectral measurements 
of higher redshift galaxies. In addition, access to spec¬ 
tral data allowed them to accurately identify galaxies 
with signihcant AGN contribution and create separate 
templates for AGN and SFGs, the latter of which is con¬ 
sidered here. The gap in spectral coverage of their tem¬ 
plates, shown as the vertical dotted lines in Figure 10, 
correspond to regions lacking spectral or photometric val¬ 
ues with which to constrain the SED, and is ignored for 
our comparison. Looking at the templates normalized 
by Ai.est(8/J.m), it is seen that the templates show re¬ 
markable agreement in SED shape for A < 24 ^m and 
lie almost entirely within the scatter in our local SED 
template. In contrast, there is clear disagreement in 
SED shape at A > 24 /rm which becomes more dras¬ 
tic at higher redshift. This is mostly due to the larger 
IRS values of these templates, which exceed the IRS val¬ 
ues observed in photometric samples at these redshifts 
(Magdis et al. 2012). For reference, the Kirkpatrick et 
al. (2012) z ^ I template has IRS = 6.5 and the z ^ 2 
template has IRS = 8.0. We associate this difference to 
the selection criteria of this sample, which required bright 
sources at 24 /im (524 > 100 /rJy) to obtain IRS spec¬ 
troscopy and which corresponds to more Air luminous 
galaxies at higher redshifts (see §6.4 for more details). 
When instead normalized by Air, slight offsets appear 
between the templates for A < 24 as a result of the 
larger Air with increasing redshift. The fact that the 
shape remains fixed, regardless of possible offsets, gives 
credibility to this technique being applicable to up to 
z ^ 2 for all bands at effective wavelengths below 24 /xm. 
In contrast, even when normalized by Air there is clear 
disagreement in SED shape at A > 24 fim which becomes 
more drastic at higher redshift. This effect, also observed 
by Magdis et al. (2012), is argued to be due to the mean 
dust temperature of galaxies increasing with redshifts. 

There are multiple physical mechanisms that give rise 
to increased dust temperatures in galaxies. Locally, sim¬ 
ilar trends are seen in galaxies with increasing values of 
Air (e.g., Rieke et al. 2009). However, for local LIRGs 
and ULIRGs there is also an associated decrease in the 
relative strength of the 8 ^m PAH feature relative to 
the FIR, corresponding to IRS values more similar to SB 
galaxies, which deviates significantly from our template. 
Instead, Magdis et al. (2012) suggest that this could be 
the result of a hardening of the radiation field, (U), in 
MS galaxies with increasing redshift {{U) oc (I -I- z)^'^®). 
Adopting Draine & Li (2007) models to fit their galaxy 
SEDs, for which {U) oc Ain/Mdust, they argue that this 
is explained by the redshift evolution of the M* — Z 
and SFR — M* relations. Another physical mechanism 
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Fig. 8.— Top: C 24 (A) conversion factor for galaxies when arranged into groups (~10 galaxies) according to L(3.6/im), a proxy for stellar 
mass, and Ljr, a proxy for the SFR. Bottom: C' 7 o(A) conversion factor for galaxies when arranged according to L(4.5/im) and Ljr. Local 
conversion factors are also shown for comparison. The dispersion in each of the groups (10 — 25%; see § 6.1) is not shown for clarity, but 
is comparable to the separation among the groups. Weak trends appear which would suggest larger conversion factors are needed for the 
higher luminosity galaxies. 


that gives rise to this effect is compactness. More com¬ 
pact star formation in galaxies can give rise to elevated 
dust temperatures and appears to occur more frequently 
in MS galaxies at higher redshifts (Elbaz et al. 2011; 
Schreiber et al. 2014). Regardless of the origin of this ef¬ 
fect, these results suggest that the 70 ^m band requires 
additional correction to be utilized as a SFR diagnostic 
as a function of redshift. We perform this analysis in the 
next section. 

Next, the templates of Elbaz et al. (2011) are consid¬ 
ered. These templates make use of redshifted photom¬ 
etry of galaxies from 0 < z < 2 to act as spectroscopic 
analog. The combination of all galaxies over this redshift 
range of results in an artificially broad FIR bump, due 
to the shifting of the FIR bump with z, and makes direct 
comparison of these templates tricky. In general there is 
good agreement in SED shape with their MS template 
and our own if this EIR broadening is taken into account. 

Lastly, the templates based on Draine & Li (2007) 
model fitting of photometric data are considered, shown 
in Figure 11. These include the templates of Magdis et 


al. (2012), and Ciesla et al. (2014). Considering first the 
z ^ 0 cases normalized by Lrest(8ttrn), we note that all 
of these model-based templates show the same excess in 
the 10 — 25 /im region compared to the spectral data 
that was seen in our own fits using Draine & Li (2007) 
models (cyan line in the Figure 11). This suggests that 
these model-based templates may not be accurately rep¬ 
resenting the intrinsic SED over this region. The region 
beyond 25 /im is likely to be more representative, as it is 
usually well fit a simple two component dust model. As 
with the spectral-based templates, the EIR bump peaks 
at shorter wavelengths with increasing redshift and also 
shows an increase in IRS. 

Taken together, it would appear that there is no 
strong evidence to suggest that the shape of the SED 
for SFGs varies significantly over the wavelength region 
of 6 — 30 /im. However, vertical offsets, corresponding 
to a constant factor offset, cannot be ruled out without 
direct comparison of SER estimates for higher redshift 
galaxies. 
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Fig. 9.— Comparison of our SFG composite template to GOALS 
photometry (U et al. 2010) and Rieke et al. (2009) templates. The 
scatter of the SSGSS template is shown as the filled gray region. 
The values of IRS = Lur/Lreat ( 8 /rm) for the Rieke et al. (2009) 
templates range from IRS = 4.8 at Ltir = IO^'^Lq to IRS = 59.6 

Ftir = lO^^X/ 0 , whereas high -2 galaxies over this luminosity 
range have IRS = 4.9 [—2.2,-1-2.9] (Elbaz et al. 2011). There¬ 
fore, there is significant FIR SED evolution that occurs for local 
(U)LIRGS that is absent at high- 2 . 

6.3. Accounting for Dust Temperature Variation 

The significant change in shape of the FIR bump makes 
the calibrations at the longer wavelengths (i.e., C 7 q(A)) 
more difficult. Comparing the SED of the higher-z galax¬ 
ies at these wavelengths to the SSGSS SED, there is a 
significant difference (up to 0.5 dex), which exceeds the 
scatter of SEDs for local SFGs. For this reason, a correc¬ 
tion to C' 7 o(A) does appear necessary if it is to be applied 
at higher redshifts. 

We correct for the dust temperature variation using the 
the SED template grids of Bethermin et al. (2012), which 
are built from the results of Magdis et al. (2012). These 
templates have been normalized by Fir. By making use 
of this grid, the observed Spitzer 70 /im luminosity as 
function of redshift is estimated while accounting for the 
changing SEDs. A demonstration of how the observed 
luminosity changes is shown in Figure 12. For exam¬ 
ple, at z = 0.5 and z = 1 the 70 /rm band measures 
rest-frame 46.7 ^m and 35.7 /rm, respectively, and the 
observed band luminosity is derived from the z = 0.5 
(purple line) and z = 1 (blue line) templates at those 
wavelengths. We perform a fit to this new conversion 
factor and present it in Table 4. The accuracy of this 
correction will be tested in the following section. 

6.4. Testing the Calibrations 

To test the utility of our calibrations, we compare 
SFRs of galaxies from other surveys to those found using 
our continuous, monochromatic values developed in this 
work. This requires a survey which has photometry avail¬ 
able in one of the calibrated bands, as well as an indepen¬ 
dent technique to measure star formation from the those 
used to calibrate our conversion factors. We choose to use 
SFRs based on Air measurements as these are the most 
readily available diagnostic for deep IR surveys. For con¬ 
sistency with our local SSGSS sample, we adopt a con¬ 


version factor of log[C'(LiR)] = 43.64 erg s“^/(M 0 yr“^), 
which corresponded to a r ^ 500 Myr constant star for¬ 
mation (see § 3.5). Furthermore, sources with significant 
AGN components need to be to identified and removed. 
It is worth noting that adopting different Ljr conversion 
factors for this analysis will only lead a constant offset 
between these two SFRs at all redshifts and that we are 
most interested in assessing where breaks from a constant 
relation develop. 

First we use the sample of 70 sources identified as SFGs 
from Kirkpatrick et al. (2012), corresponding to AGN 
contribution of less than 20%. These galaxies cover a 
redshift range of 0.3 < z < 2.5 and have full Spitzer 
and Herschel photometry. We use the Lir measurements 
of these galaxies from Kirkpatrick et al. (2012) (private 
communication), determined from IRS measurements for 
the MIR and by fitting two modified blackbodies for 
the FIR. A comparison of SFR(C' 24 (A)) to SFR(Lir) is 
shown in Figure 13, both as a function of redshift and 
Lir. There is general agreement between the values up 
to redshifts of about z ^ 1 , which corresponds to galax¬ 
ies with log[LiR/L 0 ] < 12. Given that the sources of 
Kirkpatrick et al. (2012) were required to be very bright 
in the IR to obtain IRS spectral measurements at these 
redshifts, it is likely that their sources at z > 1 are 
slightly biased to larger Lir luminosities (demonstrated 
by their larger IRS values). These values deviate signifi¬ 
cantly from deeper photometric surveys of SFGs at these 
redshifts, which is shown in Figure 13, by the templates 
of Magdis et al. (2012) (dashed cyan line). We remind 
the reader that the Magdis et al. (2012) templates are 
based on Draine & Li (2007) models, which was found to 
show significant offsets compared to the observed spec¬ 
tra of the SSGSS galaxies, and is only shown for refer¬ 
ence. We also examine the comparison of SFR(C 7 o(A)) to 
SFR(Lir), shown in Figure 14. The redshift-dependent 
correction of (770 (A) seems to work well for galaxies of 
z < 1 . 2 , for which data is available for this band. 

As a second test we use the sample from Elbaz et al. 
(2011). This sample covers a redshift range of 0.03 < 
z < 2.85 and has Spitzer and Herschel photometry. The 
Lir values for these galaxies have been determined by 
Schreiber et al. (2014), and are estimated from the Ghary 
& Elbaz (2001) template that provides the best fit to 
the Herschel data. For our analysis we only consider 
sources for which at least one photometric band covers 
wavelengths greater than 30 /im, as these cases achieve 
better accuracy of the FIR region. The photometric red¬ 
shifts of these sources are obtained from Pannella et al. 
(2014) (using the EAZY code; Brammer et al. 2008), 
and we require the sources to have suitable quality flags. 
These photometric redshifts achieve a relative accuracy 
(Az = (zphot--Jspec)/(l+-2spec)) of 3%, with less than 3% 
of cases suffering from catastrophic failures (Az > 0 . 2 ; 
Pannella et al. 2014). These sources lack spectroscopic 
measurements to identify AGN or starburst sources and 
we rely on the color-cut outlined by Kirkpatrick et al. 
(2012) for AGN and also remove sources with IR 8>8 
from our sample, which are believed to be predominately 
starburst galaxies (Elbaz et al. 2011). We follow the 
method of Elbaz et al. (2011) to determine rest-frame 
8 ^m from /c-correcting the 8 /im band (z < 0.5), 16 /im 
band (0.5 < z < 1.5), and 24 /im band (1.5 < z < 2.5) 
assuming these galaxies follow the IR SED of M82. These 
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Fig. 10. — The top of each panel shows the comparison of the composite SED of SFGs in our sample to those at higher redshifts for 
which spectroscopic information is available. The scatter of the SSGSS template is shown as the filled gray region. The SEDs have been 
normalized by Lrest(8/im) and I/jr. The sections between the vertical dotted green and red lines, corresponding to the 2 ~ 1 and z 2 
templates from Kirkpatrick et al. (2012), lack spectral data. The template of Elbaz et al. (2011) uses redshifted photometry to act as a 
spectroscopic analog. The bottom of each panel shows the residuals between our template and the other templates. The shape of the SED 
remains unchanged with redshift for A 20 /i.m, with only constant offsets occurring depending on the normalization. For A > 20 /im, 
significant SED evolution is present with increasing redshift. 
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Fig. 11.— The top of each panel shows the comparison of the composite SED of SFGs in the SSGSS sample to those at higher redshifts 
for which Draine &: Li (2007) models have been used to fit the available photometry. The scatter of the SSGSS template is shown as the 
filled gray region. The SEDs have been normalized by Lrest(8/im) and Lir. Note that when normalized by Lrest(8^m) these model-based 
template have an excess in the 10 — 25 /im region compared to the spectral data. This trend is seen in our own fits of Draine Sz Li (2007) 
models (cyan line), and indicates a limitation in the simple 3-component model typically adopted (see § 3.2). The bottom of each panel 
shows the residuals between our template and the other templates. 
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Fig. 12.— The SED of SFGs changes with redshift, as demon¬ 
strated by the templates of Bethermin et al. (2012). Taking this 
into account, the value of C 7 o(A) would change significantly from 
those derived from a. z 0 SED (dashed black line). This change is 
demonstrated as the red dashed line, which shows the 70 ^m filter 
convolution when accounting for the SED variations with redshift. 
We stress that the red dashed line is a z-dependent interpretation of 
the expected 70 ^m emission and cannot be considered a spectrum 
for an individual galaxy. 

constraints leaves us with a SFG sample of 825 sources 
with 24 /im observations and 66 with 70 fim observations. 

Using the measured Lir values from these sources, the 
same comparison is made as before and is shown in Fig¬ 
ure 15. In this case there appears to be more agreement 
among the diagnostics out to redshifts of z < 2 , with a 
Icr dispersion of 0.16 dex (45%). Small changes appear 
to develop beyond z > 2 , which is consistent with the 
observed trend of IRS going 4 to 6 by z = 3 (a differ¬ 
ence of ^0.2 dex). This trend is apparent in the ratios 
of observed 24 /rm luminosity to Lir in the templates 
of Magdis et al. (2012) (dashed cyan line). Similar to 
the Kirkpatrick et al. (2012) sample, the largest discrep¬ 
ancies occur in galaxies with log[LiR/L 0 ] > 12. Next, 
we examine the comparison of SFR(C' 7 o(A) to SFR(Tir) 
shown in Figure 16. As before, the redshift-dependent 
correction of C' 7 o(A) seems to work well for galaxies of 
z < 1, with a Icr dispersion of 0.18 dex (50%). For z > 1, 
significant differences appear but this is likely due to the 
poor sensitivity of the 70 ^m band detecting only the 
most luminous galaxies in these bands at high redshifts. 
With the limited number of sources available at z > 1, 
the reliability of our corrections cannot be determined 
for this range. 

Common to the C'24 (A)—Lir comparisons for the two 
samples considered is the trend that as one goes to z > 2 
and/or Lir > IO^^Lq the SFRs predict from Lir will 
be larger than those inferred from the MIR. This may 
demonstrate that (724 (A) is unsuitable when considering 
galaxies with log[LiR/L0] > 12 or z > 2 , but it could 
also indicate a change occurs in LiRConversion factor at 
these luminosities/redshifts. Unfortunately, it is unclear 
whether what we are observing is a redshift effect or a 
luminosity effect, as there is a degeneracy between these 
variables that cannot be resolved with our current data. 
In other words, we may simply be seeing a selection ef¬ 
fect. Given the higher sensitivity of JWST, would in¬ 


cluding sources with Lir < 12 at z > 2 follow the same 
trend of increasing IR 8 when considering the entire pop¬ 
ulation? Prior to the JWST mission, a technique that 
does not rely on Lir to determine SFRs for higher red¬ 
shift galaxies will be necessary to state confidently what 
effect is occurring. 

7. GONGLUSIONS 

We have presented continuous, monochromatic star 
formation rate (SFR) indicators over the mid-infrared 
wavelength range of 6 — 70 /rm, using a sample of 58 
star forming galaxies in the SSGSS at z < 0.2. The 
continuous wavelength coverage granted with this sam¬ 
ple has allowed for the calibration of Spitzer, WISE, and 
JWST bands as SFR diagnostics covering, continuously, 
redshifts from 0 < z < 3. We find that these diagnostics 
are consistent with monochromatic calibrations of SFGs 
in the local universe, and achieve accuracies of 30% or 
better. They also appear consistent with templates of 
high-z SFGs, with no significant evidence of variations 
in the shape of the SED over 6 — 30 /rm region. Subtle 
changes of IR 8 = LiR/Li.est( 8 / 7 m) with redshift appear 
to cause variations at z > 2 , but currently it is unclear 
whether this could be due to a selection bias at these red¬ 
shifts. Due to the signihcant changes in the FIR region 
beyond 30 fj,m with redshift, the use of this region as a 
SFR diagnostic requires correction to our local template, 
however, this has been demonstrated to work well up to 
redshifts of at least z ^ 1 . 

These powerful diagnostics are critical for future stud¬ 
ies of galaxy evolution and allows for much easier appli¬ 
cation to large survey programs with a limited number 
of MIR wavelength bands. This technique is only valid 
for SFGs, and therefore methods are required to remove 
AGN and starburst from any sample before use. With 
the upcoming JWST mission, we hope that these diag¬ 
nostics will provide important contributions as we begin 
to examine more typical main-sequence galaxies up to 
z ^ 3. 

AGKNOWLEDGMENTS 

The authors thank the referee whose suggestions 
helped to clarify and improve the content of this work. 
We also thank the MIRI instrument team for provid¬ 
ing the MIRI filter curves prior to their publication. 
AJB thanks K. Grasha, A. Kirkpatrick, A. Pope, and D. 
Marchesini for comments and discussion that improved 
the content of this paper. AJB and DC gratefully ac¬ 
knowledge partial support from the NASA ROSES As¬ 
trophysics Data Analysis Program, under program num¬ 
ber NNX13AF19G. 

This research has made use of the NASA/IPAC In¬ 
frared Science Archive, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol¬ 
ogy, under contract with the National Aeronautics and 
Space Administration (NASA). 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol¬ 
ogy under a contract with NASA. 

This work has made use of SDSS data. Funding for the 
SDSS and SDSS-II has been provided by the Alfred P. 
Sloan Foundation, the Participating Institutions, the Na¬ 
tional Science Foundation, the US Department of Energy, 









Continuous Mid-Infrared Star Formation Rate Indicators 


21 



tr 

LL 

w 


o 

cc 


CO 

O) 

o 


0.6 

0.4 

0.2 

0.0 


- 0.2 


-0.4 


- 0.6 

- 0.8 


- ★SSGSS 
" A GH-Kirpatrick-l-12 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 


0.0 0.5 1.0 1.5 2.0 2.5 0 10 20 

Zspec Number 


9 10 11 12 13 

log[L,R/Lo] 
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Fig. 14.— Comparison of SFRs estimated from C 7 o(A) and Lir for the GOODS-Herschel sample from Kirkpatrick et al. (2012). Left: 
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Fig. 15. — Comparison of SFRs estimated from C 24 (A) and Lir for the GOODS-Herschel sample from Elbaz et al. (2011). Left: 
Comparison as a function of redshift. The values show agreement for z < 2, beyond which the dataset is biased towards galaxies with 
log[Lijr/L 0 ] > 12. The distribution of the Elbaz et al. (2011) sources, along with the parameters of a best-fit Gaussian to this distribution, 
is also shown. Right: Comparison as a function of Lir. The values show reasonable agreement for cases with log[LiR/LQ] < 12. For cases 
with log[LiR/LQ] > 12, which dominate z > 2 for this sample, the monochromatic SFR is lower than the Ljr SFR. 
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Fig. 16.— Comparison of SFRs estimated from C' 7 o(A) and Ljr for the GOODS-Herschel sample from Elbaz et al. (2011). Left: 
Comparison as a function of redshift when using our calibration based on a 2 : = 0 template (filled triangles) and the ^-dependent template 
(open squares). The calibration derived from the ^-dependent template appears to work better than the z = 0 template and shows 
agreement for 2 : < 1, beyond which data is lacking. The distributions of the Elbaz et al. (2011) sources when using the 2 : = 0 template 
(filled gray) and 2 :-dependent template (open red), along with the parameters of a best-fit Gaussian to these distributions, are also shown. 
Right: Comparison as a function of Lir. 
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